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A novel method is introduced in order to treat the dissipative dynamics of quantum systems 
interacting with a bath of classical degrees of freedom. The method is based upon an extension 
of the Nose-Hoover chain (constant temperature) dynamics to quantum-classical systems. Both 
adiabatic and nonadiabatic numerical calculations on the relaxation dynamics of the spin-boson 
model show that the quantum-classical Nose-Hoover chain dynamics represents the thermal noise 
of the bath in an accurate and simple way. Numerical comparisons, both with the constant energy 
calculation and with the quantum-classical Brownian motion treatment of the bath, show that the 
quantum-classical Nose-Hoover Chain dynamics can be used to introduce dissipation in the evolution 
of a quantum subsystem even with just one degree of freedom for the bath. The algorithm can be 
computationally advantageous in modeling, within computer simulation, the dynamics of a quantum 
. subsystem interacting with complex molecular environments. 

(N ; 

\ One of the most natural ways to make a quantum system follow a dissipative dynamics is achieved by putting it 
. into contact with a thermal bath. Since usually one is not interested in the detailed time evolution of the bath degrees 
of freedom, it may also be convenient to approximate the bath dynamics by means of a classical description. When 
one faces with the problem of calculating the influence of an environment over a quantum subsystem, this approach 
leads to the representation of (a certain class of) open quantum systems 1] by means of mixed quantum-classical 
, theories. Examples can be found in many phenomena connected to quantum optics Q and quantum information 
theory Q. Typically, this is the case of cosmology where, due to the perduring lack of a full quantum theory of 
£SJ \ gravitation, one is forced to approximate formalisms in order to treat the interaction of quantum and classical degrees 
|" ■ of freedom [4]. In many situations, condensed-matter quantum systems at finite temperature can also be treated with 
mixed quantum-classical theories. In light of the above discussion, one can certainly conclude that mixed quantum- 
classical approximations can be used in open quantum systems to describe many processes which are relevant to 
various fields of research. It is worth noting that mixed quantum-classical theories [5( applied to condensed-matter 
systems can treat classical molecular bath which can be as complex as state-of-the-art molecular dynamics simulation 
techniques permit nowadays. 

t In the original constant-energy (NVE) formulation, mixed quantum-classical algorithms require many environmental 
degrees of freedom in order to describe the dissipative dynamics of the quantum subsystem. This has been shown 
within a path-integral influence functional approach 6]. However, quantum-classical dynamics has been recently 
generalized Ml in order to be unified with the constant-temperature simulation method originally developed by Nose 
and Hoover j8| (more generally, the author in Ref. Q proposed a scheme in order to unify quantum-classical dynamics 
Q-i with many energy-preserving phase space flows Q). Therefore, one can think of using quantum-classical Nose-Hoover 
(NH) dynamics in order to describe dissipative effects and, in particular, the constant-temperature relaxation dynamics 
of a relevant quantum subsystem. In practice, in order to overcome possible problems with ergodicity in classical 
phase space, it is more convenient to generalize the Nose-Hoover chain (NHC) method of Martyna and coworkers 
to the quantum-classical case and to adopt it in place of the NH dynamics. However, the choice of the NHC dynamics 
can be viewed as a mere technical point with no deep conceptual implication as far as quantum-classical theories 
are concerned. In this Communication, by simulating the relaxation dynamics of the spin-boson model [ll[, I show 
that the quantum-classical NHC dynamics can be adopted in order to describe dissipative effects in quantum-classical 
systems by means of a minimal (with respect to the number of degrees of freedom explicitly taken into account) 
representation of the classical bath. 

The quantum-classical Hamiltonian of the spin-boson model reads 

H sh = -?ina x + ( ^777 + Im^R? - CjRj&^j (1) 

where 2h~Q is the energy gap of the isolated two-state system, a z and a x are Pauli's matrices, Rj and Pj are the 
coordinates and momenta, respectively, of Nf, harmonic oscillators with mass Mj and frequencies ojj making up the 
classical bath. The other parameters of the system, i.e., (Mj,ujj,Cj), can be fixed by requiring that the harmonic 
bath is described by an Ohmic spectral density. In order to study the relaxation dynamics of this model [l2]. one 
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FIG. 1: Adiabatic dynamics of the spin-boson model: /3 = 0.3, fl = 1/3, w max = 3, £ = 0.007. The continuous black line 
represents the NVE results with Nt = 200; the white circles represent the BD results with JV b = 1 and f = 1; the black circles 
represent the NHC results with Nb = 1. 



can assume that the system is initially in an uncorrelated state with the quantum subsystem in state |up) and the 
classical harmonic bath in thermal equilibrium. The corresponding quantum-classical density matrix can be found 
starting from the full quantum expression by means of a partial Wigner transform [l3| and was explicitly written down 
in Ref. [13]. The NVE quantum-classical dynamics of this system is formally exact [l2[ (i.e., the quantum-classical 
equations of motion have the same form that would arise within a full quantum treatment) and numerical results, 
which agree very well with those obtained by means of more sophisticated path- integral iterative techniques @ , are 
available in the literature [12] . Such NVE results, which were obtained with Nb = 200, will be compared here with 
those obtained with calculations performed by means of the quantum-classical NHC dynamics and with the bath 
made up by just one harmonic oscillator (Nf, = 1). 

For the spin-boson model, the quantum-classical NHC dynamics can be defined upon introducing an extended 
Hamiltonian with a chain of just two thermostat variables 
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where T is the temperature of the bath thermalizing the quantum spin, 771, 772, p ni , p^ are the Nose variables, and 
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are fictitious masses. Following Ref. [Tj, the quantum-classical NHC bracket can be defined as: 
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where A*- NHC ' is a bracket operator whose action between two quantum-classical variables is defined as: 



(4) 
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Adopting as a convention for the point of the extended Nose phase space X = (R, 771, 772, P,Prn >Pr) 2 ), the antisymmetric 
NHC matrix reads: 
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It is worth noting that, since the Nose coordinates are intrinsically classical, a quantum-classical treatment of such 
a constant-temperature dynamics is conceptually correct and, moreover, allows one to address nonadiabatic effects. 
Other approaches [3], which do not use a quantum-classical bracket, do not seem to permit nonadiabatic calculations 
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FIG. 2: Adiabatic dynamics of the spin-boson model: /3 = 3, fi = 1/3, oj max = 3, £ = 0.1. The continuous black line represents 
the NVE results with JV;, = 200; the white circles represent the BD results with Nt = 1 and £ = 1; the black circles represent 
the NHC results with N b = 1. 



in a straightforward manner. The equations of motion can be written in the adiabatic basis as 

^X aa \x,t)=^l H Zx PP 'W) (6) 



dt 

where 
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and uj aa i — (E a (R) — E a /(R))/h. The quantum transitions operator, J, is defined as in the constant-energy case 0]. 
One wants to calculate the time-dependent quantum-classical average: 

(a, {X, t)> = J dXpfro?*' (X, t) (9) 

aa' 

where a" a (X,t) is given by Eq. ©. Details of the numerical algorithm for calculating Eq. © both in the adiabatic 
and nonadiabatic limit, can be found in Ref. [l2T |. It is useful to recall that the nonadiabatic quantum-classical 
dynamics can be pictured as a piece-wise deterministic propagation of the classical phase space point X over the 
energy surface (aa') interspersed by stochastic quantum transitions (realized by the action of J). Note also that, 
in the NVE case, one must have Nj, > 200, as proven by a cumulant expansion analysis of the influence functional 
entering the path-integral iterative procedure of Ref. Q . 

In principle, dissipation might also be described by means of the quantum-classical Brownian dynamics (BD) 
that was introduced in Ref. [15j . In such a case, the quantum-classical average of a z is still given by an equation 
similar to ^ where, however, the time evolution is achieved by means of a quantum-classical Langevin-Liouville 
operator, whose explicit expression in the adiabatic basis is known [l5T |. Such a stochastic operator is defined in 
terms of a friction constant, and of a Gaussian white noise process, £(i), with the properties (£(£)) = 0, and 
= 2ksT(S(t — t'). Therefore, it is interesting to check whether the Brownian dynamics of a bath with 
Nt, = 1 can also lead to an accurate dissipative dynamics. However, as shown in the following, especially when 
considering nonadiabatic effects, the numerical results prove that the NHC quantum-classical dynamics provides a 
scheme which is much more accurate and robust than that arising from the Brownian dynamics. 

The spin-boson model has been simulated by using dimensionless coordinates [12j j . Hundred thousand trajectories 
were produced in order to sample the initial condition in the nonadiabatic calculations [l2l |. The system parameters 
were = 1/3 and u> max = 3, while the Kondo parameter and the reduced temperature took the two sets of values 
(£ = 0.007, (3 = 0.3) and (£ = 0.1, = 3). The results of the NVE calculation, with a bath composed of N b = 200 
oscillators, were compared with those obtained with just one oscillator (iV& = 1) in the two cases where either the 
NHC dynamics or the Brownian dynamics is used. The outcome is that in the NHC dynamics a single oscillator 
provides a good numerical representation of the dissipative quantum dynamics of the spin. Moreover, it turns out 
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FIG. 3: Non-adiabatic dynamics of the spin-boson model including up to 6 quantum transitions: f3 = 0.3, SI = 1/3, w max = 3, 
£ = 0.007. The continuous black line represents the NVE result with Nt = 200; the black circles represent the NHC results 
with Nb = 1; the white circles (joined with a dashed line in order to help the eye) represent the BD results with Nb = 1 and 

C = l- 



that when nonadiabatic transitions are taken into account, the quantum-classical NHC dynamics provides very good 
results while the Brownian dynamics fails badly. Note that the the figures in this letter display the results of BD 
calculations performed only with £ = 1 (in dimensionless units). However, various other calculations were performed, 
with effective values of £ between and 10, without obtaining any improvement in the nonadiabatic case. Figures [T] 
and[2]show the results of the adiabatic calculations for the set of parameters (£ = 0.007,/3 = 0.3) and (£ = 0.1,/3 = 3), 
respectively. In the adiabatic case, both the NHC and Brownian dynamics describe well the dissipative evolution of 
the quantum subsystem interacting with a single oscillator. The inclusion of nonadiabatic transitions (up to six for 
each trajectory in the quantum-classical ensemble) shows that, with Nb = 1, the NHC dynamics is still very accurate 
while the BD evolution becomes numerically unstable at short times. This is not completely unexpected since when 
the system can switch from one potential surface to another, because of the nonadiabatic transitions, the BD dynamics 
(in the case of Nb = 1) lacks of any equilibrating mechanism. Instead, the quantum-classical NHC dynamics still 
conserves the Hamiltonian along the trajectory. Such a conservation provides a robust stabilization mechanism even 
for calculations with baths with very few degrees of freedom. Figure [3] shows the nonadiabatic results for the set of 
parameters (£ = 0.007, j3 — 0.3) while those obtained with the set (£ = 0.1, p = 3) are displayed in Fig. 0J In general, 
the results obtained with the NHC nonadiabatic evolution appear to be numerically more stable and smoother than 
those obtained with NVE dynamics. This is even more apparent in a slightly stronger coupling regime (see Fig. [4]). 

It must be remarked that surface-hopping calculations within nonadiabatic quantum-classical dynamics are, for the 
moment, limited to relatively short times because of numerical instabilities 12]. To address this issue, a quantum- 
classical non-linear formalism has been recently proposed [l6l ]. However, such long-time integration problems are 
not related to the NHC dynamics but challenge quantum-classical approximations of quantum dynamics on a more 
general level. In order to clarify this point, Fig. \5\ displays the results of a long-time calculation, performed in the 
adiabatic approximation. Since there is a great interest in the phenomenon of driven quantum tunneling [l7| , a static 
perturbation of the form H cxt — —Tvy s u z was added to the unperturbed spin-boson Hamiltonian in Eq. (fT]), and 
simulations were carried out both in the NVE (Nb = 200) and NHC case (Nb = 1) with j s /h = (l/3)fi while the 
other system parameters took the same values as in the calculations whose results are illustrated in Figs. [1] and [3] 
The results displayed in Fig. [5] shows clearly that, in the adiabatic approximation, the numerical agreement between 
the NVE dynamics (Nb = 200) and the NHC dynamics (Nb = 1) is very good even over long time intervals. 

The results provided in this Communication suggest the possibility of representing the environmental noise, leading 
to dissipative quantum dynamics, by means of deterministic NHC quantum-classical dynamics. Such an idea can 
have deep conceptual implications since there is a subtle connection between thermal and quantum fluctuations [l8[. 
Moreover, the algorithms presented here might open novel advantageous routes for the computer simulation of quantum 
dynamics in open molecular environments. Within condensed matter systems, an example that might be studied by 
means of the approach illustrated in this Communication is provided by the system recently investigated in Ref. [19( : A 
retinal chromophore molecule evolving according to short-time quantum coherent dynamics in bacteriorhodopsin. As 
already shown in Ref. [2(| , a minimal model of such chromophore-protein systems can be built by explicitly considering 
the chromophore molecule itself and the nearest-neighbor amino-acids, belonging to the tight-binding pocket in which 
the chromophore is contained. On the short-time scale of the coherent quantum dynamics of the chromophore, 
one might think of representing the dissipation entailed by the rest of the protein by means of a deterministic NHC 
dynamics with a minimal bath. It is also worth to note that there is a great interest in the field of quantum information 
in the phenomenon of driven quantum tunneling j 1 71 ] - In particular, recent work focuses on time-dependent external 
driving [2l|. In order to deal with such situations, one should generalize the quantum-classical approach presented 
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FIG. 4: Nonadiabatic dynamics of the spin-boson model including up to 6 quantum transitions: /3 — 3, Q — 1/3, o; max = 3, 
£ = 0.1. The continuous black line represents the NVE results with Nb = 200; the black circles represent the NHC results with 
Nt, = 1; the white circles (joined with a dashed line in order to help the eye) represents the BD with Nb — 1 and ( — 1. 




FIG. 5: Adiabatic dynamics of the spin-boson model under the action of a static perturbation of the form H cx t = —^ a (J z : 
P = 3, fl — 1/3, oj max = 3, f = 0.1, 7s/ft = (l/3)f2. The continuous black line represents the NVE results with Nj, = 200; the 
black circles represent the NHC results with Nb = 1. 



here in order to unify it with the techniques of non-equilibrium molecular dynamics simulations. Although non-trivial 
algorithmic issues might be expected, this is something which is possible in principle and that deserves a thorough 
future investigation. 

In conclusion, the quantum-classical NHC dynamics may find interesting applications in various fields: In fact, it 
may be used to simulate not only systems in the fields of chemical-physics or biophysics but also in quantum-optics 
and quantum computing [3j. 
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